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We present a Fermi liquid model for the overdoped and optimally doped cuprate supercon- 
ductors. For the normal state, we provide an analytic demonstration, backed by self-consistent 
£ — , Baym-Kadanoff (BK) numerical calculations, of the linear in temperature resistivity and linear in 

1/energy optical conductivity, provided the interacting Fermi liquid has strong peaks in its density 
of states (van-Hove singularities in 2 dimensions) near th e chemical potentia l fi. Recent ARPES 
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I. Introduction 



expts. by Valla et al., Science 285, 2110 (1999), and e-print cond-mat/0003407 directly support the 
D ' linearity of the one-particle scattering rate everywhere in the Brillouin zone hereto obtained. We 

»' | , show that the origin of this linearity is the linear in energy term of the imaginary part of the carrier 

■ susceptibility. Moreover, we verify that the interactions tend to pin the van-Hove singularities close 

to (i. We show that the low energy dependence of the susceptibility can have a purely fermionic ori- 
gin. We introduce an ansatz for the susceptibility of the carriers, which we postulate to be enhanced 
in an additive manner due to the weak antiferromagnetic order of the CuC>2 planes. Inter alia, this 
ansatz may explain the appearance of the spin resonance peak (observed in neutron scattering) in 
■""^ . the normal state of the cuprates. Further, we obtain particularly high transition temperatures T c 

from our BK-Eliashberg scheme by using this ansatz: we have a d x 2_ y 2 gap with T c > 120°K for 
O \ nearest neighbour hopping t — 250meV. 

PACS numbers: 74.72.-h, 74.20.-z 
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The nature of the many-body state of the cuprate superconductors is a central question for the understanding of 
these materialsErtlflE.g. one long standing puzzle has been the elucidation of the origin of the linear in temperature 
in-plane resistivity!] and linear in 1 /energy optical conductivity observed in the optimal doping regime and to a good 
extent in the overdoped regime. Vice-versa, the answer to this question should shed light on the character of the 
carriers and, subsequently, on the-.superconducting transition. Here we address these issues based on a minimum 
unconventional Fermi liquid modelEl. Our model comprises strong peaks in its density of states (van-Hove singularities 
in 2 dimensions) near the chemical potential. We show that it accounts in a natural, comprehensive and internally 
consistent manner for several normal state characteristics. The introduction of an ansatz for the susceptibility of the 
carriers further allows us both to propose an explanation for the origin of the spin resonance peak and to obtain 
particularly high <i-wave transition temperatures T c . Overall, our results make a strong case for a Fermi liquid 
approach to the optimally doped and overdoped cuprates. 

We perform a combination of analytical and numerical many-body calculations in the context of our model. The 
rest of this paper is organized as follows. In Section II we write our many-body approximation, which we used in 
our numerical calculations. We emphasize both that our treatment is relevant for overdoped and optimally doped 
cuprates (see Section VI for the underdoped regime) and that the results presented in Sections III-V depend only 
quantitavely and not qualitatively on the specific Hamiltonian and approximation thereof etc. In Section III we discuss 
our analytical and numerical results for the linear in max(T, e) scattering rate of the carriers, in connection with the 
existence of van Hove singularities (vHs) close to the chemical potential. We also show that a conductivity linear in 
(1/T, 1/e) follows. In Section IV we discuss both the fermionic origin of the energy dependence of the Millis-Monien- 
Pines susceptibility and our ansatz of eq. (^6|) for the susceptibility of the carriers. We show that this ansatz may 
explain the appearance of the so called spin resonance peak, seen in neutron scattering experiments, in the normal 
state of the cuprates. In Section V we discuss the superconducting transition in the frame of the ansatz. Finally 
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Section VI contains a summary of our results. In the Appendix we examine the role of doping- induced disorder on 
the carrier susceptibility. 



II. General framework 



We assume that we deal with a Fermi liquid, albeit an unconventional one, as will become apparent from our 
discussion of the scattering rate of the carriers below. We choose the 2-dimensional Hubbard Hamiltonian as a 
specific model for our numerical calculations (c.f. the last paragraph of Section I on this). 

H = e k cl,a C k,a + 2]y2 C i+ q ,A>- q ,-* C k\-<?Ck,a • (1) 

cj, is an electron creation operator and is the electronic tight-binding dispersion suggested by angle-resolved 
photoemission (ARPES) experiments - e.g. see! - and LDA calculations! 

£fe = — 2£(cos k x + cos k y ) — At' cos k x cos k y — 2i"(cos 2k x + cos 2k y ) . (2) 

It is assumed here that the lattice constant is equal to unity and is the same along the two crystal axis a, b in the 
planes - hence k x ,k y e [—%, n]. N x N is the discretization of the Brillouin zone. n 

We consider the fluctuation-exchange diagrammatic approximation (FLEX) of Bickers, Scalapino and WhitcEl for 
the Hamiltoniap, which consists in summing bubble and ladder diagrams. FLEX is a Baym-Kadanoff conserving 
approximatiorB, meaning that there is a free energy functional $[G] of the Green' function G, such that the self- 
energy £ is given by the relation £ = S$[G]/5G. We thus obtain a set of self-consistent equations for G(k,e n ) and 
S(fc,e n ) : 

G(k, e„) _1 = G c (fc, e^) -1 — E(fc, e„) , (3) 

Go(fc,e„) = ^ , (4) 

ie n + (i- e k 

T,(k,e„) = -—^^V{q,uj m )G(k-q,e n -uj m ) . (5) 
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The potential V(q,uj m ) is given by 



V(q,u m ) = V ex (q, uj m ) - V H {q,u m ) , (6) 

^^' r^Tvi ■ (8) 

The susceptibility Xo{Q: w m) is given by 

Xo(q,u m ) = ~(T/N 2 )J2G(k + q,e n +Lu m )G(k 7 e n ) . (9) 

£„,!; 

fi is the chemical potential and the Matsubara frequencies are e n = (2n + 1)ttT and u) m — 2mirT for fermions and 

bosons, respectively. We solve numerically this self-consistent set of equations, working with a given number M of 
Matsubara frequencies and discretization of the Brillouin zone (M = 256 — 480 and N > 64). 
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There have b§en ajrumber of similar numerical calculations on the normal-phase and the superconducting transition 
of the cupratesa'tjiU, with FLEX being a particularly popular approach. 

All the convolution operations ace done by using the Fast Fourier Transform (FFT) , in order to cut down calculation 
time. We use Pade approximantsE3 to analytically continue our results to the real frequency axis. 



III. On the scattering rate of the cuprates 



We have analytically obtained a scattering rate linear in the maximum of the temperature T or energy e, for a 
Fermi liquid with strong density of states peaks - van-Hove singularities (vHs) in 2-d - located at an energy e v n close 
to the chemical potential fi. 

The derivation relies on the relation (||) for the self-energy, which is valid quite generally in the frame of a BK 
approximation, irrespectively of±he specific Hamiltonian and approximation thereof. See also the discussion following 
eq. (po|). It can easily be showntj that ImT,(k, e) is given by the following formula at finite temperature : 

im£ fl (fc,e) = ^2lmG R {q,e- Lo)ImV R {k - q,u) {coth(cj/2T) + tanh((e - w)/2T)} . (10) 

Taking 



we obtain 



ImG R (k, e) = -nS(sk,e), Sk,e = e + £t - e k - ReT,(k, e) , (11) 



/m£ K (fc,e) = -tt> ImV H (k - q,s q , e ) {coth(s q , t /2T) + tanh((e - s 9j£ )/2T)} . (12) 



Setting ImG R (k,e) equal to a delta function is a reasonable approximation for this purpose, in view of the typical 
sharp spike feature of ImG R (k,e) shown in fig. 1 - also see figs. 2 and 3, all of which are representative of our 
numerical solution of eqs. (3)-(9). Further, numerically ImG R (k,e) is very small compared to the band energy for 
small couplings, and the difference of IrnS(k, e), as seen in our numerical calculation, for small and large coupling 
constants is mostly quantitative rather than qualitative. 
We write 



V q {2n+1) (0) X 2n +1 

^ (2n + l)! 



ImV R (q, x)=^ 1 ^y iM , (13) 

n=0 



where V q n \o) is the n— th derivative of ImV R (q, w = 0) with respect to w. This is true for an electronically mediated 
interaction, with a polarization which is a regular function of ui (see also eq. ( p0| ) below) . There are only odd powers 
of bj in the series because the imaginary part of the susceptibility is an odd function of energy - e.g. prf. eq. (2.63) of 
Pines and Noziereai3. One possible exception to this is given by Gonzalez, Guinea and Vozmediandl3. The authors 
showed that for the underdosed LSCO-type Fermi surface (FS) and for momenta K connecting two inflection points 
of the FS, the imaginary part of the susceptibility goes like |o; | 1 / 4 in 2 dimensions. However, this fact will influence 
the final result for the scattering rate only for a small range of momenta k satisfying k = K + q a - c.f. eqs. (14) and 
(15), and q is given in the paragraph below. Moreover, we have already emphasized that our picture is not valid in 
the underdoped regime. 

First we consider the low T limit. The sum over q is dominated by the van- Hove singularities at the points q . 
Assuming that c v h — e <? + < ReT,(q a , e) > [this relation is misprinted in the journal version of the paper] is close to 
/i, the tanh has a vanishing contribution at the vicinity of e q + ReT,(q, e) ~ fx (note that for e q + ReT,(q, e) < /i and 
e q + ReYi{q 1 e) > fi + e the contributions of tanh and coth annihilate each other in the low T limit). Hence 



ImZ R (k, e),-,EE " ,o , .If , (14) 



For sufficiently small Vq n \o)^ Vn > 1, we obtain 



(2n + l)! 
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a* = * 2^ ^-,(°) » * (2n+l)! ' ( 5) 

where c = fi — e v u- This relation is valid for e + c < e c , where the latter is the characteristic energy beyond which the 
infinite sum on the right becomes comparable to the Vq term. Also, for energies beyond the bandwidth W (W = 8t 
for the non-interacting system), Imx , and hence ImV (see below), decay to zero. These considerations yield the two 
energy crossovers 

ei = \fi-e v it\ , £2 = min{e c + e V H - V, W + e v n — m} j ( 16 ) 

while the assumption above for a k leads to 

ImY. R (k,e) ~ -a k {e + c) , e x < e < e 2 ■ (17) 

For e > e 2 im£ gradually decreases, due to the finite bandwidth of the system. Finally, we note that if the Fermi 
surface approaches a van-Hove singularity at q , v£_ (0) should become bigger, being pronortional to l/(ye kF kp) (as 
implied by the standard Fermi liquid result for the imaginary part of the susceptibilityEO - see the discussion below 
on the susceptibility of the cuprates). _ . 

We consider now the high temperature limit T > (fx — € v h)/4&3. We see immediately that 

ImE R (k,e) = -nY,ImV R (k-q, Sq , t ){2T/s q ^ + 0{ Sq ^/2T)} . (18) 

(Note that the term of order T of this sum is reminiscent of the left-hand side of the sum rule - c.f. Pines and Nozieres0 
- limq^o Jo du)Imx (q,u)\e(q,u})\ 2 /u> = —Nir /mc 2 s , with N being the total particle number, c s the speed of sound, 
m the effective mass, and e(q,u) the dielectric function.) The sum is dominated by the van-Hove singularities at the 
points q Q , thus yielding 



ImS fi (fc, e) ~ -2Ttt ^ ^(0) = -2a k T . (19) 
q~q, 



Here we made use of the condition above for a k - In addition, it is straightforward to see from our analytic treatment 
that Im'S R oc x 2 , x =max{T, e}, when both T, e — > 0. In all respects we have a genuin e Fermi liquid. 



Note added: ARPES expts. by Valla et al., Science 285, 2110 (1999), and preprint |cond-mat/0003407| , have very 
recently shown that in optimally doped P^Si^CaC^Os+y the one-particle scattering rate is linear in max{T, e} over 
most of the Fermi surface, in support of our picture. „ 

A brief comment here. It has been known long ago - see e.g.El - that the scattering rate becomes linear in T for 
T > ll>b /4, with lvb being the characteristic boson frequency mediating the carrier interaction. Our treatment shows 
that u>b here is nothing else but the fermionic energy \x — e v h- 

The prefactors in the r.h.s. of eqs. (fD]) and ( |l9| ) differ by a factor of 2. This is in agreement with experiments in 
YBa2Cu307_5 and Bi2Sr2CaCu208+j;Ej^ where the factor is found to be in the range 2.1 - 2.2 . We note-that the 
"marginal Fermi liquid" phenomenology of Varma, Littlewood, Schmitt-Rink, Abrahams and Ruckensteir£3 gives a 
factor of 7r instead. 

We emphasize that the T and e dependence of the result are independent of k - thus leading necessarily to a 
linear in T resistivity and a linear in 1/e optical conductivity, even with inclusion of vertex corrections in the 
calculation. The reason for this being that {T, e} are obtained as overall prefactors, for the relevant T and e regimes, in 
such calculations. E.g. the Kubo formula yields a(w) = (e 2 /u) J2 k e v l G(k, e + uj)G(k, e) [1 + S(k,e)] [f(e+u) — f(e)], 
where v k is the group velocity, S includes vertex corrections from the Ward identity, and / is the Fermi occupation 
factor. The main T and uj dependence in the integrand is in the one-particle self-energy in G and in /. Doing the k 
sum, we get the dominant contribution from the poles of G. Now, the one-particle scattering rate is linear in max(T, e) 
everywhere in the Brillouin zone, and this linear dependence appears in the denominator of a. Here we assumed that 
the vertex corrections do not have-.a strong temperature dependence over a substantial part of the Brillouin zone. 
Indeed, Kontani, Kanki and Uedao have recently shown numerically, in the frame of the FLEX approximation, that 
vertex corrections ate, small for the resistivity, and do not change its T dependence. 

Hlubina and RiceEil considered analytically a model of interacting fermions with a vHs close to [i. However, they 
find a scattering rate similar to ours only close to the vH region, and different otherwise. As a result, their resistivity 
goes like T 2 ln 2 (l/T). In their 'hot' and 'cold' spots scenario, relying on strong scattering off antifcrromagnetic 
fluctuations, they obtain an average scattering rate similar to ours, and numerically a linear in T resistivity (however, 
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they seem to assume that the group velocity is finite along the whole Fermi surface - cf. betweep-eqs. (2.5) and (2.6) 
ofEJ) . Similar results are also obtained in the antiferromagnetic scenario of Pines and Stojkovict3. 

A note on phonons. As they form a - presumably small- part of the effective potential V, they provide necessary 
momentum dissipation, yielding a finite resistivity. However, the linear T dependence of the latter is not specifically 
influenced by phonons in our model. 

Returning to the derivation of the scattering rate above, we observe that the overall behavior of ImV closely follows 
Imxo{q,uj), as 

ImV{q,u>) = Im Xo (q,w) Hq^)\ 2 , (20) 



Imxo is odd in cj, while |e| 2 is even. Eq. (20) follows from any screened interaction between the carriers. Hence the 
argument for the linear in energy and temperature behavior of t _1 (T, e) is equally generic. It relies essentially on a 
large coefficient for the linear in energy term of ImV - i.e. of Imxo - and the presence of van-Hove singularities near 
the Fermi surface. The result holds regardless of the dimensionality of the system. However, it is important that a 
significant part of the spectral weight be included in the strong peaks of the density of states lying close to fj,. 

What is more, in our self-consistent numerical solution we observe that the energy e v u of the singularities is pushed 
by the interactions close to the chemical potential - see fig. 4. This result is especially pronounced when we use the 
ansatz for the susceptibility of the carriers of eq. ( p6| ) below. Then we find typically for n ~ 0.87 — 0.95 and for a 
broad range of t' , t", U 

M - evH < t/20 . (21) 

The shape of the self-energy £(£;, e) of the interacting system is responsible for the modification of the density of states 
iV(e), through the relation N(e) — —Tr ImG(k, e)/ir. A trend for the transfer of the spectral weight is indicated by 
the fact that JmE(fe, e) has a peak below /i and a dip above it. T-hp numerical result concerning thc-approachment 
between e V H and \x has been known for some years. Si and LevinEl and Newns, Pattnaik and Tsueic3 observed the 
pinning of the vHs close to [i by using a U — > oo mean field r^ave boson approximation of a model with Cu 3d 
and O 2p orbitals. Recently, Gonzalez, Guinea and VozmedianoEZI were able to obtain analytically the essential part 
of the approachment between the vHs and /x with a first order renormalization group treatment in the context of 
the Hubbard, model. A review of related work in the frame of the so-called van-Hove scenario has been given by 
MarkiewiczEa. This pinning of the vHs close to fi seems to be a plausible explanation for the common characteristic 
of a good many cuprates whose van-Hove singularities are located between 10-30 meV below the Fermi surfaceEJ (see 
also the next section). 

It is interesting that the electron doped Nd2- :c Ce a ;Cu04+ ( 5 which has a van-Hove singularity much below the Fermi 
surface, i.e. at approximately /x-350 meV, as shown by ARPESEj, has a usual Fermi liquid t _1 (T) = const. T 2 ln(TjEl. 
This lends support to the picture described above. Along the same line, the resistivity of T^P^CuOetH (Tl-2201) 
switches over from linear to quadratic with increasing doping from the optimal to the overdoped regimeE3, which we 
suspect to be an indication of the vHs moving well away from fi. _ . 

Finally, in further support of the relevance of the vHs in the transport properties of the cuprates, Newns et. alL^I 
and Mcintosh and KaiserEil have shown that the thermal conductivity of the cuprates can be well accounted for if the 
vHs are located very close to the Fermi surface, as discussed above. 

The numerical solution of the many-body system always corroborates our analytical result for the self energy at 
finite temperature. im£(fc,e) turns out to be essentially linear in energy in the interval t\ < e < 62. A linear 
dependence of r lmE(fc, e) as a function of either T or e was also obtained in the numerical work of Beere and AnnettE3, 
Kontani et al.E3 and Si and Levino. Note that in figs. 2 and 3 we show the self-energy for the set of the system 
parameters which yields the highest transition temperature T c , if use of the ansatz of eq. (]2^) is made. The linearity 
of JmS(/c, e) with e is even more pronounced for other combinations of i', t" and n. iro£(fc, e) has always the correct 
parabolic Fermi liquid bevahior for e — ► 0. Furthermore, the energy interval of linear behavior expands as the energy 
e V H of the (extended) van-Hove singularities at the (vicinity of the) points q = (±7r, 0), (0, ±7r) approaches [i. 

Another feature of the density of states as seen in our treatment - cf. fig. 4 - is the following. The non-interacting 
density of states has two minor peaks at the bottom and top of the spectrum respectively (in fig. 4 the top one is 
a vHs). As the strength of the interaction increases, these two peaks are washed out, as a result of the self-energy 
which becomes substantial in magnitude for energies away from /i - cf. figs. 2 and 3. 

At the moment it is not clear whether the present mechanism of the linear scattering rate can explain the experimen- 
tally obsepasd T 3 dependence of the Hall resistivity of the cuprates. A way to explain it has been found by Stojkovic 
and PinesEj, using an electron interaction peaked at Q — (±tt, ±7r). Their argument can be slightly modified, so that 
it works for our form of the electron potential V - given by eq. (@) - butr-with a modified effective % peaked at Q, as 
we propose in the next section - cf. eq. (Eg) and below. Kontani et al.c3 have shown that vertex corrections in the 
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frame of FLEX have a drastic influence on the T dependence of the Hall resistivity, in marked contrast to the case of 
the longitudinal resistivity. 



IV. On the susceptibility of the cuprates 



The low energy dependence of the susceptibility of the cuprates. The Millis-Monien-Pines susceptibilit 

xmmp(q,u) = 2 . — ^4 — ~i — ' ( 22 ) 

1 + 4 (<? - Q) - iu/u S f 

has been used to fit the low energy part of the susceptibility of the cuprates in both NMR rate and inelastic neutron 
scattering (INS) experiments. Here Q = (±7r, ±tt). The short range antiferromagnetic (AF) order, a remnant of the 
parent antiferromagnetic materials, with correlation length £, is responsible for the peak of the susceptibility for q 
near Q. Typically £ is of the order of the lattice constant (£ decreases as the doping increases, and e.g. £ ~ 2 for 
optimally doped YE^CusOy-^), while ujsf « 10 — AOmeV. 

The origin of the small magnitude of ^s^p-has remained elusive thus far. E.g. Sachdev, Chubukov and Sokol have 
interpreted it as a damped spin wave modeEZI. Spin waves are clearly observable in underdoped cuprates. However 
to date there is no experimental proof that they are strong enough in the normal phase of the optimally doped and 
overdoped regimes. The proximity of the system to an antiferromagnetic instability, i.e. Vq xiQ^) — 1) where V q 
and u>) are some appropriate coupling and susceptibility respectively, can in principle explain the small magnitude 
of ujsf, as pointed out by Millis, Monien and PinesE3. 

Here we propose that an alternative explanation - which may coexist with the latter - is the following fermionic 
origin for ujsf- 

First iinwever, let us present a proposal for the susceptibility which has been put forward by Onufrieva and Rossat- 
MignocH. This starts by viewing the CuC>2 planes as a lattice of plaquettes centered on the copper site with four 
nearest neighbour oxygen sites. A Hamiltonian TL, reminiscent of but more comprehensive than the one of the t — J 
model, was introduced in terms of the Hubbard operators. In this formulation, the itinerant carriers which propagate 
via Cu spin flips are clearly separate objects from the localized Cu spins with short range AF order. A diagrammatic 
approach was developed in the frame of TC, leading to the following RPA-type total susceptibility 

/ n XAF{q,uj) +XF(q,u) 

Xt(q,v) = — —, -. r— -. rr . (23) 

1 + Jq{XAF(q, w) + XF(q, w)) 

J q is the effective Cu spin exchange interaction, xf is a purely fermionic susceptibility and XAF(q^) is due to the 
localized spins. Xt(g,w) encompasses in an appealing way the idea of the entangled carrier-spin dynamics in the 
cuprates. Furthermore, this approach is able to account to a good extent for the variation of the total susceptibility 
as a function of doping and temperature, as measured by INS. 
Now, we use the result for xt above with 

Xi ? 2 

XAF(q,u) = „ --2 jj-^- , XF{q,u) =XFo{l + iu/uj + 0(uj 2 )) ,q^Q. (24) 

Let us suppose that f{uj) = iuj/ujs- If ujs 3> uj q and Jq£, 2 Xi < 1; an d taking XF{q,w) = Xo{q,w) (as given by eq. 
@); we essentially recover XMMp(q,u) - which is itself an approximate form of the true susceptibility - with 

ujsf -> u{q) = ■ -^f — - . (25) 

uj + uisJqXt (1 + <r (<Z - Q) ) I U Xi) 

From the numerical solution of our system, we easily obtain values of lj comparable to the experimentally relevant 
ones, when the van-Hove energy e v n is near /x, with ui a scaling quickly towards zero as \i — — > 0. Hence ui can 
be interpreted as oj {qp) = Ve qF qp - c.f. the non-interacting Fecpii liquid result u> {q) — vpqtJ. The small difference 
A 1 — e vH is observed in a good number of cuprates. E.g. in refo there is a compilation of several cuprates, the van- 
Hove singularities of which are located between 10 - 3CLmeV below the Fermi level (c.f. the discussion in the previous 
section). Also, Blumberg, Stojkovic and Klein (BSK)B suggested that this characteristic may be true irrespective of 
the doping, as long as the latter is appropriate for superconductivity. This is based on ARPES experiments on the 
bilayer YBa2Cu307_a. ARPES remains the best diagnostic probe for the Fermi surface of the cuprates. Yet it has not 
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proved possible to perform measurements on many other compounds, especially the monolayers such as T^E^CuOg+s 
etc. The point here is the following. By fitting the ARPES data BSK show that one of the two effective bands - the 
anti-bonding one - formed by hybridization of the two layers by interlayer coupling has a chemical potential only some 
20 - 50 meV above the van-Hove singularity at (0, 7r), irrespective of the doping regime. It is then clear that these 
carriers, with a large density of states, give rise to a small uj as discussed above. Hence it is very interesting to know 
how universal this band-structure characteristic of the cuprates is, as it may explain naturally the magnitude of to. 
Furthermore, it would be interesting to determine experimentally, e.g. by INS, the value of a) for Nd2- a; Ce£ ! ;Cu04+i. 
In that case, to should be enhanced as a result of the van-Hove singularities being far away from the Fermi surface. 
In the Appendix we discuss the (non)influence of weak disorder on the value of w Q . 

The antiferromagnetic ansatz for the carrier susceptibility. We thereby propose that the effective non- 
interacting susceptibility (i.e. without interaction lines connecting the particle-hole lines) of the carriers is given by 
the following ansatz 

Xo(q,u) -^xt ff (q>u) = Xo(q,u) + axAF(q,u) . (26) 

Xo is given by eq. (9) above, xaf is the antiferromagnetic susceptibility of the localized Cu spins and < a < 1 is 
a dimensionless weighting factor, which in principle depends on doping, band structure, temperature (presumably a 
decreasing function of the latter) etc. The intuitive idea is that the carriers shoulcLperceive the ordered antiferro- 
magnetic background of the CuC>2 planes, even in the absence of phase separationLj|_J. In this way, the fermionic 
susceptibility acquires an antiferromagnetic enhancement, which may then influence the effective carrier potential and 
pairing, through an electronically mediated interaction - see also the next section. Furthermore, x% can explain the 
temperature dependence of the Hall resistivity of the cuprates, as we mentioned at the end of the previous section on 
the scattering rate. Finally, we note that with this effective xt^ our many-body approximation remains conserving, 
since the relation £ = SF/5G between the self-energy, the free energy and the Green's function is still pifid, as xaf 
does not depend on G. We note that this is consistent with the work of Onufrieva and Rossat-MignodEa mentioned 
above. The carriers and the localized spins form two distinct, albeit interrelated, systems. The additive form of the 
ansatz is also compatible in spirit with xt(q,^) of eq. (|23|). _ 

We mention here the alternative treatment of the spin and charge susceptibilities of Imada, Fujimori and Tokuracj. 

Taking the ansatz of eq. (|2^) at face value in the context of our many body scheme means that both the charge and 
spin susceptibility of the carriers acquire an AF enhancement. Currently we cannot prove this ansatz. We emphasize 
that our ansatz can be taken to apply only for the spin susceptibility of the carriers. In that case, in the frame of 
our Baym-Kadanoff scheme, only Xo(q^) entering the ladder diagrams with opposite particle-hole spins would be 
replaced by x% ■ Quantitatively, the differerence between this case and the case in which both charge and spin 
susceptibilities are enhanced is small (for relevant values of the AF enhancement) - c.f. the discussion on the critical 
temperature T c in the next section. 

In our numerical implementation, we consider two similar forms for the AF susceptibility, namely (A) 

4 

xi F {q,Um) = x ^r- 1 9(uj c - \ Um \), (27) 

i=l 

with Fj = + (q — Qi) 2 , 0J C being a cut-off, and (B) 

Br \ v n ^ m ~ (2w m /7r)arctan(w m /L>) - TiD + (2TiD/ir) arctan^) 

x A F\q^m) - a v ^2 _ (r. t £>) 2 ' ^ ' 

with Qi = (±7r, ±7r) and D being a cut-off frequency, above which ImxAF(q,u) = 0. Form (B) has appeared i 
Here the characteristic spin wave frequency obeys ll>s ok. £~ z , and the z = 2 scaling regime has been assumed, iu. 
agreement with the analysis of Sokol and Pines for the optimally doped and overdoped regime of Y-F^CusOv-a etc.cfl 

The spin resonance peak of the cuprates. INS experiments in YBa2Cu307_a - see e.g.c3 and therein - and 
Bi2Sr2CaCu208+jJl3 have revealed the existence of a strong peak in the spin triplet channel, centered at Q = (±tt, ±tt) 
and at a characteristic resonance energy uir ~ 30 — 40 meV. Although this peak is usually seen exclusively in the 
superconducting state, it has been observed up to temperatures ~ 250°K in YBa2(Cuo.gg5Zno. oos^OyEj for ujr t-W^X 
meV. Interestingly, this fact cannot be accounted for by most of the theoretical models so far available - see.e.gM3'c£l 
for refs. - as these models require the onset of superconductivity. An exception is the model of Buluto, which, 
however, differs drastically from ours. 

The use of the ansatz of eq. (^6|) in the spin channel only may account in a natural way for the appearance 
of the resonance peak in the normal state through a bilayer effect. Both YBa2Cu307_5 and Bi2Sr2CaCu20s+j, are 
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bilayer materials. Here one can define bonding and antibonding bandsO-a e±k — ek±t±(k), t±(k) being the fc-dependent 
interlayer hopping element. Further, one defines the susceptibilities Xo±(q, k>) = Tr[G + (+q, +uj)G±+G-(+q, +uj)G t ]. 
When e V H — A* is small, Xo-{q,^) has a narrow peak at q = Q and lo ~ A, ssihere A is the bilayer splitting at 
q ~ 9o = (tt, 0), (0, 7r), i.e. the van Hove neighbourhood. Liechtenstein et al.tj have shown in the frame of the 
FLEX approximation that A is drastically reduced for a finite interaction, compared to the non-interacting value 
A c = 2t^(q = q Q ), such that it becomes comparable to the experimental lor. Our ansatz can be taken to apply to 

Xo±, yielding Xo± -> xtj = Xo± + a Xaf- As a result X-(l, u) = xl- (q, <*>)/(l - Uxl-(q,u)) is strongly peaked at 
Q for an energy oj c r ~ A, and may account for the experimental observations. Of course we require Xo- > Xo+ for 
this to work, which is valid for a range of the parameters t' ,t" ,n. 

One can ask the question: why does the resonance peak not appear in the normal state in general? As demonstrated 
here, the amplitude of the AF enhancement needs pto be sufficiently large for the peak to be visible. Zn is known to 
enhance AF fluctuations in the Cu02 planes - e.go, which if interpreted as yielding a larger (a xaf) contribution 
in Xo \ is i n agreement with the results above. On the other hand, as we discuss in the next section on the 
superconducting transition, too strong a factor (a xaf) reduces T c - c.f. T c — 93°K for the pure YBa 2 Cu 3 07 versus 
T c = 87°K for the Zn doped material mentioned above (also see the last paragraph of section V). 



V. On the superconducting transition of the cuprates 



With the solution of the normal system at hand, we solve the gap (Eliashberg) equationca 

A(fc, e„) = ~ V p( k - k '> e " ~ e ») G ( fc '' e «) G (- fc '> -OA(fc', O . (29) 

k',e' n 

This form of the equation is valid close to the transition temperature T c only. The pairing potential V p is given by 

V p (q,LO m ) = V x (q,Wm) + Vh(q,(V m ) , (30) 



V ^ m)= l-U* U xl{q,. m ) ' (31) 
V h {q,. m )= ^T m \ ■ ^ 

1 - Uxo{q,u m ) 

We have assumed that the gap is an even function of both its momentum and energy arguments in order jm. write 
Vp(q,u> m ) in this form. We also assumed that the gap is spin singlet, as Knight shift-experiments have shownatl. The 
superconducting state of the cuprates is probably a generalized BCS state r^ee e.g.EEl - with the transition being due 
to a momentum anisotropic potential, as envisaged by Kohn and LuttingeiO. _ 

The symmetry of the gap is determined by the exact shape and sign of V p . Based on symmetry groundsel as well 
as on experimental evidence, we expect a d x i_ y i or a s wave gap - but see alsoLa. (In principle, higher order even 
angular momentum harmonics are also possible.) The highest T c 's correspond to a d x 2_ y 2 gap, and are obtained by 
including an antiferromagnetically enhanced susceptibility in the calculation, following our ansatz of eq. (p6|). In 
passing, let us note that in general the proximity of e V H to /i plays a less significant role than the ansatz in raising fCr 
- but c.f. below. Nevertheless we also obtain a s-wave gap under the same conditions, but with a much lower T'^l. 
This is consistent with the experimental situation. Most of the cuprates appear toriiaKe, a d x 2_ y 2 gap at T c - e.g. c.f.E. 
Experiments point to the opening of a secondary order parameter gap at T' c <C T c E-3t3. 

We note here that the van-Hove singularities at q a tend to suppress a s-wave gap if the pairing potential V p >0 
(V p defined above is negative for sufficiently large and negative U and/or an appropriate phonon coupling, and 
positive otherwise). The (likely) s-wave gap of Nd2_ 2; Ce 2; Cu04+ ( j could originate from e.g. the fact that the van-Hove 
singularities are 350meV below the Fermi surface, and hence ineffective here, or from V p <0 for some relevant parts 
of the phase space in this material, or possibly from both facts. Yet another possibility which may coexist with the 
above is that V p is strongly peaked close to zero momentum owing to the band structure of this material. All these 
factors can result in the s-wave T c being higher than the d-wave T c , and hence in the dominance of the former channel 
over the latter. 
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We obtain the following (near) optimum set of parameters for the d-wave T c . We take t = 250meV and £ = 1 
(increasing the latter leads to a reduction of T c , see below). For form (A) - eq. (^) - we obtain T c ~ 125° K, for 
al„ = 0.2eV _1 , uj c = 4i, t' = -O.lli, t" = 0.5t, f7 = 1.27eV and n = 0.88. For form (B) - eq. (§|) - we obtain 
T c ~ 105°^, for aX„ = 0.5eF -1 , I? = 32t, t' = -O.lli, t" = 0.25t, 17 = 1.47eV and n = 0.91. The variation of T c in 
i4T with £ is as follows for this last set of parameters : (94, 1.5), (91, 2U/90, 3), (88, 5). 

Monthoux and Pines obtained similarly high transition temperaturesc9 with their approach, which uses V(q, u) = 
g 2 XMMp(q,<jj) as the effective carrier-carrier potential. Their approach has no room for the Hubbard U and the sub- 
sequent screened carrier-carrier interaction though. Overall, our approach yields significantly enhanced T c 's compared 
to the standard FLEX-type treatments!!! 

We obtain an optimum value of n for the following reason. \o - and hence x% - is a decreasing function of n. We 
allow the coupling to increase up until b — U\ e J $ {Q,lo = 0) saturates to a value close to and below unity. Then V 
and V p (given by eqs. (6) and (32) respectively) are increasing with n. And so does the characteristic scattering rate 
entering G. Thus the optimum value of n for the highest T c corresponds to the overall highest kernel V p \G\ 2 in the 
gap equation (p9|). 

Further, we have done an extensive search of the parameter space to locate the optimum parameter set for the 
highest T c . From our data it appears that the variation oLX c as a function of the system parameters is smooth and 
that the optimum parameter set above is a global optimurr£3. 

Assuming that only the spin susceptibility acquires an AF enhancement according to our ansatz, yields a T c which 
is lower by 9% for £ = 1, but only lower by 2% for £ = 3, if we make use of Xaf' with the optimum parameters above. 
The variation of T c with £ is as follows here : (96,1), (92,1.5), (90,2), (89,3), (88,5). 

There is an optimum value of a for the maximum attainable T c . This is again due to the form of the pairing 
potential V p above : the AF instability condition allows Ux e J^ < 1 only. Now, for a given value of the latter product, 
the highest V p - which in principle yields the highest T c as well (c.f. above) - will correspond to the highest possible 
U. This in turn corresponds to a smaller X% i and hence to a small but finite optimum a. In the next paragraph we 
discuss relevant experimental evidence. .— . 

Zheng, Kitaoka, Ishida and AsayamaEJ made a very interesting empirical observation. Namely, among the hole 
doped cuprates, the highest T c 's correspond to a combination of both optimum total carrier concentration n x 2_ y 2 +2n PtJ 
in the planes (n x 2_ y 2 being the concentration of holes of Cu-3d orbital character and n Pa of 0-2p orbital character) 
as well as of a reduced (probably minimum) imaginary susceptibility - as deduced from NMR experiments. In fact, 
Zheng et al. noticed that the highest T c 's are obtained for a reduced ratio n x i_ y i /2n Pcr ~ 1 ; and, moreover, that such 
a trend is correlated with a reduced relaxation rate (1/Ti)c„ oc T lim^^o Y] q Imx(q, w)/w, in properly normalized 
units (here A q is the hyperfine coupling). This last fact points to a reduced Imx(q,^) etc. These conclusions are 
in accordance with our picture, which yields both special values of the filling factor n - for this also c.f. e.gflO - as 
well as special small values of the product (a xaf) as a prerequisite for the highest attainable T c 's (also c.f. the last 
paragraph of section IV). 



VI. Summary 



To summarize, we present a single plane Fermi liquid model which for the normal state can explain the salient 
transport properties, the low energy dependence of xmmp, and their relation to the existence of van-Hove singularities 
close to the Fermi surface. E.g. we obtain analytically a scattering rate linear in max(T, e), within appropriate T and 
e bounds, for all momenta in the Brillouin zone. This result yields directly a linear in T resistivity and linear in 1/e 
optical conductivity. The introduction of an ansatz for the susceptibility of the carriers allows for an understanding 
of both the appearance of the spin resonance peak in the normal state and the temperature behaviour of the Hall 
conductivity. Further, by using the ansatz we obtain significantly enhanced d x i_ y i wave transition temperatures T c . 
Attention is paid to Nd2- x Ce :E Cu04+i5, the properties of which, despite appearances, we believe to be fully consistent 
with those of the majority of cuprates. 

In brief, let us discuss the possible connection to the physics of the underdoped cuprates. Strong experimental 
evidence-suggests that they are in a phase separated regime, with AF domains of spins separated by stripes of 
holesEjEl. One can envisage that with doping increasing towards the optimal regime, the stripes melt into an 
effective Fcpni liquid, and the physics described here is recovered. Further models on the underdoped cuprates can 
be found inEl 

The author has enjoyed discussions with Yia-Chung Chang, Gordon Baym, Joseph Betouras, Girsh Blumberg, 
Antonio Castro Neto, Lance Cooper, Sasha Liechtenstein, Peter Littlewood, Jorg Schmalian, Raivo Stern, Qimiao Si 
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and Branko Stojkovic. This work was supported by the Research Board of the University of Illinois, the Office of 
Naval Research under N00014-90-J-1267 and NSF under DMR-89-20538. 



Appendix - Disorder effects in the susceptibility of the carriers 

Returning to the origin of a small u> , which was discussed in section IV, another option would in principle be the 
disorder inherent in the cuprates. The dopants are randomly positioned in the crystal structure, thereby creating 
an effective disorder potential for the carriers in the planes. We calculated the effect of disorder by considering 
the dopants as isotropic point scatterers, with a density rij=l% and typical scattering strength V s = — 1 eV (i.e. 
< 8t). These parameters give a residual impurity scattering rate less than 2 meV , consistent with experiments on the 
cuprates. For the calculation of the susceptibility we used the diffusorH. We only used a band structure with t" = 0, 
and only took into account non-magnetic disorder. The non-interacting Green's function now becomes 

G' (k,e n ) = — 1 , (33) 

ie n + (j, - e k + a l {e n ) 

with Oi(e n ) = riiV s / (1 — (V s /N 2 ) G (k, e„)). The susceptibility is given by 

^ LI - riiV/P{q,uJ m ;e n ) i 

and P(q, uj m ; e n ) = (1/N 2 ) G(k + q, e n + Lo m )G(k, e n ). However we found no evidence, for the parameters above, 
of uj being influenced by disorder. 
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FIGURE CAPTIONS 

Figure 1. Full Green's function G R (kp,e) for (a) kp = ^ h ° undar y along (tt,0) — > (tt,tt) and (b) kp — k'^ la 9 onal along 
(0, 0) -> (tt, tt), for t = 250meV, t' = -O.llt, t" = 0.25i, U = 1.5eV, n = 0.91, at T = 105°^. 

Figure 2. Self-energy Y.(kp,e) for kp = j t b ° undar y along (w, 0) — > (tt,tt) (continuous line) and kp = f t <^ a 9° nal along 
(0, 0) — > (tt, tt) (dashed line), for the same parameters as in fig. 1. (a) depicts liaS R (kp, e) and (b) depicts KeT,(kp , e). 
A (quasi)linear energy dependence of ImS R (kp, e) can be seen for energies below 0.5eV - also c.f. text. 



11 



Figure 3. Self-energy T,(kp,e) for kp = ] i b ° undar y along (ir,0) — > (7r,7r) (continuous line) and kp = k F la9ona along 
(0,0) -> (7r,7r) (dashed line), for t = 250meV : t' = -O.llt, t" = 0.25i, U = lA7eV, n = 0.91, at T = 105°K. The 
carrier susceptibility has an antiferromagnetic enhancement according to the ansatz of eq. ( p6| ) here (see section IV) , 
with £ = 1, aX Q — O.beV^ 1 and D — 32t. This is the optimum T c case when using form (B) - eq. H2a ) - of our ansatz. 
(a) depicts ImS fl (fcf,e) and (b) depicts KeT,(kp,e). A (quasi)linear energy dependence of ImE fl (fcjr,e) can be seen 
for energies below Q.5eV. 

Figure 4. Evolution of the density of states. Dashed line : non-interacting system, with parameters as in fig. 3. 
Continuous line : same system with U = 0.8ey. Notice the transfer of the central van-Hove peak towards the chemical 
potential, the disappearance of the two satellite peaks - see text, and the broadening of the total spectral width. 
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